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The statistics of photons emitted by single multilevel systems is investigated with emphasis on 
the nonrenewal characteristics of the photon-arrival times. We consider the correlation between 
consecutive interphoton times and present closed form expressions for the corresponding multiple 
moment analysis. Based on the moments a memory measure is proposed which provides an easy way 
of gaging the non-renewal statistics. Monte-Carlo simulations demonstrate that the experimental 
verification of non-renewal statistics is feasible. 
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The arrival times of photons emitted by single quan- 
tum systems have become a task of routine measurements 
[lH9|. Several methods are currently in use for the anal- 
ysis of the recorded photon time traces, such as the sec- 
ond order field-correlation function [1-y, LLO, [ll| , photon- 
number statistics |12| . Il3| ] , exclusive and nonexclusive in- 
terphoton probability density functions (PDF) [3l. ll4 - |lrj| . 
or Mandel's Q-function [3, II H E3 ■ Usually, the PDF 
of photon-arrival times is considered to depend solely on 
the arrival time of the previous photon, assuming tacitly 
that photon emission is a renewal (semi-Markovian) pro- 



cess [1J| . Consequently, multiple interphoton time PDFs 
are factorized and cast into products of one-interphoton 
time PDFs [3|, [l4|. In contrast, a nonrenewal process 
indicates a memory, since the photon-arrival time PDF 
depends not only on the arrival time of the previous but 
also on the arrival time of the photon before last, and 
consequently on the particular realization of the previ- 
ous photons' time trace [la, LUa, ll 



For ensembles of microscopic photon sources the pho- 
ton statistics is expected to be renewal, however, ex- 
periments reported on single and coupled quantum dots 
Is [l9rl21| , on single pairs of coupled molecules [8 L 22J, |23| , 
and on two-state dynamics of single molecules [7j could 
be considered for the investigation of nonrenewal proper- 
ties in the photon statistics. Recently, a renewal indicator 
was introduced for the study of conformational fluctua- 
tions of single molecules [24J . This indicator is closely re- 
lated to Mandel's Q-function and relies on the statistics 
of the number of photons recorded in a given time inter- 
val. In this paper we consider another technique which is 
based on the correlation between consecutive interphoton 
times. We apply a multiple-moment analysis of consec- 
utive interphoton times and propose a measure M for 
deviations from renewal statistics. 

The time evolution of a multilevel quantum system 
interacting with the radiation vacuum can be given in 



terms of the reduced density matrix p by (h = 1) |ll[|l5| 

P = Cp = i[p, h] + J2 m (sr P s+ - \ [s+s-, P ] + ) ,(1) 

where the Liouvillian C consists of the Hamiltonian H , 
which includes the interaction with the classical driving 
field, and of dissipation in the Lindblad form. As usual, 
S^~ ($£ ) are lowering (rising) operators for the i-th tran- 
sition, "fij denote for i = j the spontaneous emission 
rates and for i ^ j cooperative decay rates deviating 
from zero if the difference of the two involved transition 
frequencies is smaller than the inverse radiation-bath cor- 
relation time, \ui— u>j\ < 1/t c [18]. Assuming the rotating 
wave approximation, C does not depend on time so that 
the evolution of p is given by p(t) = e Ct p(0), where p(0) 
denotes the state at time zero. 

For the description of state collapses upon photon 
detection, several theoretical approaches, pioneered by 
the Monte-Carlo wave function technique [25|, were de- 
veloped. These approaches rely on quasi-continuous 
photon-emission measurements to introduce system 
states conditioned on whether a photon is detected or 
not |18L 1261) . Accordingly, the Liouvillian is split into two 
terms 1181 



C = C c + K 



(2) 



where C c governs the time evolution of the conditioned 



and non-normalized density matrix p c (t) 



Xot 



P(0) 



U c (t)p(0), subject to a zero-photon outcome of the mea- 
surement up to time t. The second term in Eq. ([2]) rep- 
resents the collapse (reset, recycling) operator 1Z [lSJ] to 
reset the density matrix upon a photon detection event 
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(3) 



The dimensionless detection efficiency rj is introduced to 
account for the fact that a state collapse takes exclu- 
sively place when the emitted photon is also detected 
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FIG. 1: Multi-level systems under consideration: (I) three- 
level A-system, (II) cascade three-level system, (III) four-level 
system motivated by a pair of interacting two-level systems, 
and (IV) a two level system jumping stochastically and ra- 
diationless between two states. Heavy arrows for laser-light 
driven and spontaneous radiation transitions and light arrows 
for radiationless transitions. 



16]. According to Eq. ([2]), also the conditioned Liouvil- 
lian C c = C — 1Z depends on r\ in a unique way. 1Z op- 
erating on p c (t) at random times generates a stochastic 
process and thus the average survival probability Pq (t) 
of no-photon detection up to time t is [11| 



P (t) = Tr{p c {t)} = Tr{U c (t)p } 



(4) 



provided that a photon was recorded at time zero. Cor- 
respondingly, po is the average state just after photon 
detection and is given by normalizing the collapsed sta- 
tionary state, po = 1Zp ss /Tr {lZp ss }, where the stationary 
state satisfies Cp ss = 0. 

Recording the times of state collapses generated by re- 
peated application of the operator 7ZU c (t) mimics the 
time traces of photon detection in a particular single 
quantum system experiment. The conditional density 
matrix right after the n-th photon detection of a sequence 
of exclusive detection times {ti, • ■ ■ , t n } with ti > t%—\ is 
then given by 



p c (*i,*2, ..,*»)= [r+Yl^Kudu 



fa_i 



Pa 



(5) 



where the time ordering operator 7+ ensures that the 
operator at the latest time is on the far left. The trace 
of p c {ti,t 2 , ..,t n ) in Eq. ([5]) provides the detection PDF 
of a particular time sequence 

p n (Ti,T 2 ,---,T n ) = Tr{\T + Y[ i=i KUc(Ti) paj , (6) 

where Ti = ti — tj_i are intcrphoton times. Furthermore, 
the PDF P2(t) of detecting a second photon at time t, 
given a detection event at any previous instance, follows 
from summing up all possible realizations of two consec- 
utive interphoton times [151 ] 



ft(*) 



P2(t-Ti,Ti)dri 



(7) 



Generally, referring to Eq. © the PDF P n (t) of the 
n-th photon at time t results from the n — 1 fold con- 
volution of the operator TUA c (t), where for completeness, 
Pi(t) = Pi(t). For the quantitative analysis of p n , we 
examine the moments to order mi , i = 1 , ■ • • , n for n con- 
secutive detection intervals. These moments can readily 



be calculated using a moment generating function tech- 
nique in several dimensions 

P mi ,..,m n = ( II*=l / dTl ^ ) P™( Tl ' ' • ' ' T «) 

= Tr{(r+n" =1 (-l) (mi+1) ^^C (mi+1) ) Po} ,(8) 

where, recalling, the ordering operator 7+ ensures that 
the operator at the latest time is on the far left. Eq. 
(jHJ allows for an easy numerical calculation of multi- 
ple moments. In case of renewal, TZp(i) does not de- 
pend on t, in other words the state after a collapse is 
independent of the state just before the collapse. Con- 
sequently, Pn(Ti, ■ ■ -T n ) of Eq. ([6]) can be factorized in 
terms of the one-interphoton time PDF, P^(t\, • • • , r„) = 
nr=iPi( T j)> where the superscript R denotes renewal. 
Furthermore, the arrival PDF of the second photon is 
P^(t) = Pi(i) *pi(t), where * indicates convolution and 
multiple moments reduce to products of individual mo- 
ments 



rmi ,-",7717, 
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(9) 



Differences between the PDFs p n (Ti, ■ ■ ■ ,r„) and 
Pn(ri, ■ ■ ■ , T n ), Pi{t) and P^(t), or between the moments 
of Eqs. (jHJ and J9]) may be used to demonstrate whether 
the initial state po is recovered after photon emission and 
the process is renewal or whether the state-resetting de- 
pends on the current state and the process is nonrenewal. 
The deviation from renewal is a signature of the lack of 
information about the system state after photon emission 
and indicates the memory present in the correlation be- 
tween consecutive photon-arrival times. Envisaging the 
experimental verification of NRS we concentrate on two 
consecutive time intervals and propose the following mea- 
sure 



M = pi.\l\x\ - 1 



(10) 



which can be determined directly from the experimental 
time traces and can easily be predicted using Eq. flSJ). M 
takes on both signs and an analysis of bi- valued waiting- 
time sequences indicate that tentatively M is negative 
when shorter and longer waiting times are likely to oc- 
cur alternatingly and is positive when both, shorter and 
longer waiting times are likely to be bunched. 

The multiple moments of Eq. ((8} and the measure 
M. of Eq. (fT0|) represent the main result of this paper 
and we consider A4 as a characteristic quantity comple- 
mentary to other statistical measures, e.g. the photon 
coincidence probability (PCP) g^(0). We also studied 
the covariance of the functions P^it) and P% (t), however, 
such an analysis requires binning of the time traces and 
is therefore less direct for gaining information about the 
memory in the photon statistics. 

To illustrate the NRS in photon counting we consider 
a representative set of level schemes, shown in Fig. 1. 




FIG. 2: The three-level cascade system, a) Memory M, PCP 
<7 (0), and fluorescence-excitation intensity Ip as a function 
of the laser detuning A. The arrow indicates the A value 
used in b) and c). b) Renewal and nonrenewal second- photon 
arrival PDF P^{t) and P2(t), respectively, and one-photon 
arrival PDF Pi(t). c) Two-photon arrival PDF P2(i~i, T2). d) 
Memory M as a function of the detuning A and the Rabi 
frequency Q. 



For the three-level A-system (I) renewal applies, because 
once a photon is detected the system is reset to the same 
mixed state, no matter the state before emission. In this 
respect, the A-system equals a two-level system which 
always collapses to the ground state. In contrast, NRS 
arises for systems (II-IV): For the three level cascade (II), 
upon the detection of a spectrally unresolved photon, the 
populations and coherences of the states |1) and |2) be- 
come proportional to the ones of |2) and |3) prior to emis- 
sion, respectively [18j . In the four level cascade (III), the 
reset operator ladders populations and coherences down 
the levels: |4) -> (|3), |2)) -» |1). The non-cascade sys- 
tem (IV) shows a TLS flipping stochastically between 
two states. The flipping may be associated with changes 
of spectral and dynamical properties so that bunching of 
short and long interphoton times and thus NRS results. 
Summarizing, if the state after emission depends on the 
state prior to emission, the photon time traces obey NRS. 
The PCP can be discussed accordingly, namely, for sys- 
tem (I) it is zero, and is non-zero for systems (II) and 
(III). However, for system (IV) the PCP is zero although 
the process is nonrenewal in general. 

We next discuss the level schemes (II-IV) in more de- 
tail. The Hamiltonian in Eq. ([1]) can be written as 

H = ffo + E* h^( S t + 5 0, wher e fi* is the Rabi fre- 
quency of the i-th transition resulting from the interac- 
tion with the classical driving field. For level scheme (II) 
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FIG. 3: Jumping two-level system, a) If scaled as indicated 
and M as a function of the laser detuning, similarly as in Fig. 
2a. The arrow indicates the detuning used in Fig. 2d. b) A4 
as function of v and K, on linear and log scales, respectively, c) 
M as a function of the detection efficiency r\ on log-log scales 
for two values of fC. The dashed line for the r]~ 2 dependence, 
d) Second-photon arrival time PDFs P^(t) and Pi(t). e) M 
and confidence intervals M ± om as a function of the number 
of detected photons. Wiggly lines for Monte-Carlo simula- 
tions and dashed lines for predictions. The dash-dotted line 
gives the predicted value M — 1.04. A = v in panels b)-e). 



we write for the zero-order Hamiltonian 



H, 



(ii) _ 



= -6iS?Sr-(6 1 + 6 2 )SZSi 



(11) 



where 57" = |1)(2| and S 1 ^ = |2)(3| are lowering operators 
with rising operators defined accordingly. S\ = Wl — W21 
and 62 = ujl — ^32 are differences between the laser fre- 
quency wl and transition frequencies Wij. Results are 
shown in Fig. 2 for the parameters: i?i = Q% = fi = 
27; 7n = 722 = 7)712 = 721 = 0,w 32 - W21 = 87, and 
77 = 1. In Fig. 2a the excitation fluorescence inten- 
sity I F = Tr {Kp ss }, the PCP 5 ( 2 '(0) = Tr {ft 2 // 8 } /If, 
and A4 are compared as functions of the detuning A = 
h(Si + 5 2 ). The fluorescence shows a maximum located 
approximately at resonance with the the lower transition 
(<5i ~ 0), followed by a peak at A ~ 0, where coherent 
two-photon absorption and cascade emission are likely to 
occur. The intensity at resonance with the upper transi- 
tion ($2 — 0) is weak because of weak pumping of level 
2. The PCP indicates photon antibunching in the range 
of the lower transition and photon bunching in the range 
of two-photon absorption and of the upper transition. In 
agreement with the above discussion, A4 takes on nega- 
tive values in the range of A ~ 0, where alternating short 
and long interphoton times are probable. At resonance 
with the upper transition, A4 is weakly positive indi- 



eating minor bunching of short and long waiting times. 
P%(t) and P%{t) displayed in Fig. 2b deviate consider- 
ably from each other and similarly a strong asymmetry 
is apparent in the PDF P2( T i7 T 2) upon interchanging t\ 
and Ti demonstrated in Fig. 2c. In Fig. 2d M. is mon- 
itored as a function of A and ft. A minimum close to 
A ~ and Q ~ 27 is clearly visible. We have found 
that \M\ drops roughly as rj~ 2 so that the experimental 
verification of NRS requires a high detection efficiency. 

Motivated by recent investigations of pairs of identical 
and interacting quantum systems [8J, [2l|, [23 1 we studied 
system (III). Depending on the parameters the results 
(not shown here) were similar to those of system (II) 
which is obvious when states 2 and 3 are superposition 
Dicke states [2fi [28( , so that the behavior is governed by 
the 1-0-2 and 2f->4 transitions or lo3 and 3-H-4 transi- 
tions. 

We finally report on the non-cascade, jumping two- 
level system (IV) where the ± states indicate for instance 
two different molecular, or lattice nuclear configurations 
[7j, |l7[ , spectral diffusion of ultracold molecules in con- 
densed media [29], or two different spin configurations. 
The dynamics is described by extending the density oper- 
ator, p = (p_,p + ) T , and correspondingly the Liouvillian 
and reset operators [IT | 



from the condition 
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(12) 



where /C± denote the jumping rates between the two 
states and where C± and 1Z± are the Liouvillian and 
resetting operators of the two states. For illustration we 
assume the full symmetric case where only the transi- 
tion frequencies are different for the two sates. Thus the 
system is described by the Hamiltonian H± = (— A ± 
v)S± S± + hQ(S± + S±), where A and v are laser de- 
tuning and frequency displacements from the transition 
center, respectively. Furthermore, j± = 7, 1Z± = 1Z, 
and JC± = /C. Numerical results are shown in Fig. 3 
for the parameters Q = 2j,v — 27,77 = 1,}C = 7/IOO, 
except when they appear as variables or are specially in- 
dicated. Ai is positive throughout all calculations and 
peaks close to the resonances of the two states. Fig. 3b 
indicates large positive M. for K. <C 7 and for |A| ~ \v\. 
Fig. 3c shows how the 77 dependence of M. crosses over 
to the asymptotic r\~ 2 behavior and how the crossover 
is shifted to lower values of r\ with decreasing JC. The 
second-photon arrival PDFs P^(t) and P%{t) in Fig. 3d 
differ only at longer times which indicates that these 
quantities are not appropriate for providing evidence of 
NRS. 

To demonstrate the experimental feasibility of measur- 
ing M, we report Monte-Carlo simulation results llll25|. 
Choosing a random number r distributed uniformly in 
[0,1], the detection time t n of the n-th photon follows 



Ti{U c {t n -t n -i)p c (ti,---,t n -i)} = r 



(13) 



where p c is the normalized conditioned density matrix of 
Eq. (5). By resetting and normalizing the state at t n , 
the initial state of the next interphoton cycle is obtained. 
Fig. 3e shows, how M. converges as a function of the 
photon number to the predicted value. Also presented 
are confidence intervals M ± <jm , which are estimated 
from the variance, assuming statistical independence of 
the moments: a 2 M = a 2 t i /fif + 4(/ui,i//if ) 2 er^ x , where 

a-% = M2 - Mi and a 2 ^^ = p 2 ,2 - Mi,i- 

The experimental investigation of the NRS requires 
the measurement of two consecutive intervals, so that 
the arrival times of three consecutive photons have to 
be recorded. Depending on the time scale, this can be 
achieved using a single detector, however, for time scales 
shorter than the detectors' dead time, at least three de- 
tectors are needed. A comprehensive description of ex- 
perimental data has to account for the detectors' dead 
time and for the ubiquitous background photons. 

In conclusion, we have shown that NRS is plausible 
in the fluorescence of multi-level systems and that the 
indicator A4, proposed for the identification of the non- 
renewal property, is experimentally feasible. For cas- 
cade systems M. may be small at low detection efficiency, 
so that advanced experimental techniques are required, 
while for non-cascade multi-level systems M. may be 
large also at low detection efficiency. 
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